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High-mass X-ray binary systems are powered by the stellar wind of their donor stars. The X-ray 
state of Cygnus X-1 is correlated with the properties of the wind which defines the environment of 
mass accretion. Qiandra-HETGS observations close to orbital phase allow for an analysis of the 
photoionzed stellar wind at high resolution, but because of the strong variability due to soft X-ray 
absorption dips, simultaneous multi-satellite observations are required to track and understand 
the continuum, too. Besides an earlier joint Chandra and RXTE observation, we present first 
results from a recent campaign which represents the best broad-band spectrum of CygX-1 ever 
achieved: On 2008 April 18/19 we observed this source with XMM-Newton, Chandra, Suzaku, 
RXTE, INTEGRAL, Swift, and AGILE in X- and 7-rays, as well as with VLA in the radio. 
After superior conjunction of the black hole, we detect soft X-ray absorption dips likely due to 
clumps in the focused wind covering > 95 % of the X-ray source, with column densities likely to 
be of several 10^^ cm^^, which also affect photon energies above 20keV via Compton scattering. 
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Introduction 

One of the currently still puzzling issues in our understanding of microquasars is the exis- 
tence of the different states (low/hard, high/soft, very high/steep power-law, . . . ), which distinguish 
themselves by the soft X-ray luminosity, the X-ray spectrum, the timing properties, and also the 
nature of the jet and its radio spectrum. One assumes that these states with their in many ways 
different properties correspond to different configurations of the accretion flow which cause dif- 
ferent channels or efficiencies of energy dissipation. For a high-mass X-ray binary system, the 
global framework of mass accretion onto the black hole is defined by the powering wind of the 
donor star. One possibility is that changes in the wind properties (like density, velocity, ionization 
parameter, . . . ) might finally trigger state transitions. 

This contribution focuses on the wind of the well-known persistently bright X-ray source 



CygnusX-1. Its companion star HDE226868, a ~18M0 09.7 type supergiant [|11|], fills 97% 
of its Roche lobe [^]. The stellar wind is therefore strongly focused due to tidal interactions ^. 
The detailed structure of this wind defines the environment of mass accretion onto the black hole on 
the largest scale. At a luminosity of Lo.5-250keV = 3 x 10^' erg s^^ (which is obtained already in the 
low/hard state |^]), the X-ray source produces a considerable feedback on the wind by photoion- 
ization [Q], which can create a complex wind structure within the binary system. The parameters 
inferred from the Ha emission line are (anti-)correlated with the X-ray flux: in the low/hard state, 
the wind is probably denser and therefore less strongly ionized, which allows for a more efficient 
wind acceleration and thus a lower accretion rate onto the black hole |^. 

We have obtained multi-satellite observations of CygX-1 in the hard state close to superior 
conjunction of the black hole, i.e., at phase 0orb ~ in the 5.6 day binary orbit ||3|, |], which are 
ideal to investigate the wind: at the low inclination / k, 35° no eclipses of CygX-1 by its 
companion star can be observed. Nevertheless, absorption dips, which are likely related to the 
focused stream covering the X-ray source, are often detected close to <^oxh = [|l|. 

1. The Chandra and RXTE observation of 2003 April 

1.1 Advantages of a multi-satellite observation 

In the first part of this contribution, we show some of the results from a simultaneous obser- 
vation of CygX-1 performed by Chandra and RXTE in 2003 April shortly before 0orb = 0. (More 
details are given by Hanke et. al [Q, 10].) The light curves shown in Fig. [T] clearly show the ad- 



vantages of a multi-satellite observation: Due to RXTE's low Earth orbit, there are data gaps of 
~45 minutes length in each RXTE-orbit (lasting ~90 minutes) from passages through the South 
Atlantic Anomaly (SAA). The high variability which is present on timescales of hours and which 
is related to the stellar wind could hardly be investigated with the RXTE data alone. The (continu- 
ous) soft X-ray light curve taken by Chandra is shaped by absorption dips. Only the joint spectral 
analysis of the non-dip continuum allows to combine the strength of the different instruments. For 
example, the PCA spectrum clearly requires an Fe Ka emission line at 6.4 keV, but the resolution 
is not good enough to decide whether it is narrow or relativistically broadened. In the simultaneous 
Chandra spectrum, a narrow component is detected (see Fig. but its flux is not sufficient to ex- 
plain the PCA feature, from which we infer the presence of a broad iron line as well. (More details 
on iron lines in Cyg X-1 are given by M. A. Nowak in his contribution to this volume [jl^.) 
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Figure 1: Light curves of CygX-1 during the simuhaneous Chandra / RXTE observation of 2003 
April. Top: i?XT£-HEXTE (20-250keV). Middle: RXTE-PCA (4-20keV). Bottom: Chandra -HETGS 
(0.5-12keV order events). The dashed lines define the non-dip times. 



1.2 The non-dip spectrum 

The non-dip spectrum can be obtained by filtering for high count rates, see the dashed horizon- 
tal line in Fig. [I[ The results of this subsection are thoroughly described by Hanke et al. (2008) [^]. 

The Chandra -HETGS spectrum of the bright source CygX-1 shows neutral absorption sig- 
natures of the interstellar medium (and probably also source environment). We detect structured 
absorption edges, namely the Fe L2 and L3 edges at 17.2 and 17.5 A (0.721 and 0.708 keV) - 
which J. C. Lee describes in detail in her contribution in this volume -, the O K-edge at 22.8 A 
(0.544 keV) with the Ka resonance absorption line, and the Ne K-edge at 14.3 A (0.867 keV) with 
the Kj3 resonance absorption line. These complex absorption edges can be modelled with tbnew\ 
an improved version of tbvarabs, which contains those fine structures. 

The non-dip spectrum is dominated by absorption lines from the stellar wind of highly ion- 
ized ions, namely H- and He-like O, Ne, Na, Mg, Al, Si, S, Ar, Ca, and Fe, as well as many Fe 
L-shell absorption lines, see Fig. ^. We fit these features with Gaussian absorption lines and infer 
equivalent widths and Doppler shifts with respect to the rest wavelengths of the identified transi- 
tions. In a second approach, we have defined a fit-function (in ISIS [|l^, 15]) for the whole series 
of absorption lines of an ion. Theoretical Voigt profiles were used which take into account the 
oscillator strengths of the strongest transitions. The ion's column density is a direct fit-parameter 
of the model, avoiding an explicit curve of growth analysis. The velocity shifts inferred from both 
methods are rather low, mostly |v| < 200km s^^ The latter result can be understood within a sim- 
ple model for the spherical wind, e.g., a generahzed CAK-model [p], [l4|] : 



v(r) = v(r) • - with v(r) = vq + (v„o-vo)- 1 

r \ r 



(1.1) 



See http : //pulsar . sternwarte . uni-erlangen . de/wilms/ research/ tbabs. 
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Figure 2: Qiandra-HETGS CygX-1 non-dip spectrum with identification and fits of absorption lines of 
H- and He-like ions (Lyman and Helium series; black labels) and Fe L-shell ions (blue labels). Note the 
narrow and broad Fe Ka line (red) at 1.94 A (=6.4 keV, see second x- (energy-) axis on top of each panel). 
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Figure 3: Wind velocity, projected against the 
black hole. Absorption Unes with low velocity 
Vjad can be observed from the dark purple re- 
gion close to the black Thales circle of Vrad — 0, 
v'lad > can only be observed from the red re- 
gion inside this circle, and Vrad < only from 
the blue region outside of it. The brighter the 
red or blue, the higher is |vrad('^)|. The dashed 
gray lines are contours of constant Vjad in steps 
of 200 km s"^ The highlighted area with the ^o^b 
scale shows the region in the orbital plane with 
I cos I < sin /, which equivalently corresponds 
to all possible lines of sight (for / — 35°). 



For the velocity observed in X-ray absorption lines, one has to consider the projection against 
the black hole, Vrad('^) = cos a{f) ■ v(r) , where a{f) is the angle between wind velocity v(r) and 
the line of sight towards the black hole. Figure ^ shows Vrad (f) as a colorscale plot. On the cir- 
cle passing through the center of the star and the black hole, a = 90° by Thales' theorem, and 
thus Vrad = - independent of the assumed velocity law, as long as it is spherically symmetric. 
The circle separates (possible) absorbing regions of the wind seen at a redshift (inside) from those 
seen at a blueshift (outside). The location of absorbers seen at a low velocity is thus well con- 
strained. It is important to realize that this result does not depend on the model for the wind 
velocity, provided that v(r) points radially away from the star. 

Due to the system's inclination, not every line of sight in the orbital plane can actually be 
observed: the angle t? between the inclined line of sight / = (sin/, 0, cos/) and the binary axis 
a = (cos(27r0orb); sin(27r0orb)) O) is given by cos ^ = 1 o a = sin/ • cos(27r^orb)- The range of 
possible t> values is obviously limited by a low inclination. (For a face -on geometry with / = 0°, 
there is clearly no orbital variation, i.e. t? = 90°, while every < t? < 180° is possible for edge-on 
systems with / = 90°.) The region with | cos < sin / is highlighted in Fig. || between the crossing 
solid gray lines. For a spherically symmetrical wind, this region in the orbital plane is equivalent to 
all lines of sight experienced by an inclined observer. The green labels show the position of these 
lines of sight for given orbital phases 0orb> where (^orb and 1 — <porb are identical. 

As a result, the observation of low Doppler velocities in absorption lines in front of the X-ray 
source constrains the location of the absorber - rather independently of the wind model. This 
location is in good agreement with XSTAR simulations of the photoionization zone. The number 



density of H-atoms, «h> is obtained from the velocity profile (Eq. 1.1) and the continuity equation, 
M^^ = uniu ■n}i{r) -Anr^ ' Here, pL = 1.4 is the mean molecular weight per H-atom and 



niH is the mass of an H-atom. Inserting the numbers Voo = 2100 kms and 7?^ = 11 Rq [11] gives: 

nn{r) = l.lxlO^^ cmT^ ■ d{r /R^) , where d{x) = ^ (1.2) 

voK + (l-vo/voc)-(l-l/x)'^ 

Between the stellar surface and the black hole, x = r/R^, varies between 1 and 2.4 and the average 
wind density is ~ 10^' cm^^ for vo/vc» = Vtherm/vt» = 0.01 and a reasonable j3 f« 1. 
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Figure 4: The Mg xi ttiplet of 
resonance (r) absorption and in- 
tercombination (i) as well as for- 
bidden (f) emission lines is de- 
tected with two pairs of i- and 
f-lines in the non-dip Cban- 
dra-HETGS spectra of Cyg X- 1 . 
The red-shifted pair stems from 
an emitting plasma component 
of higher density and can be 
identified with the focused wind. 
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While these previous findings from the non-dip spectrum could be explained by a spheri- 
cal wind only, there are also signatures of a focused wind: Figure ^ shows the Mgxi triplet 
of Is^(^So) ls2p(^Pi) "resonance" (r), ls2(^So) ^ ls2p(¥i,2) "intercombination" (i) and 
Is^(^So) <— ls2s(^Si) "forbidden" (f) transitions. The i- and f-emission lines show up as two 
components each: one pair is almost at rest - in accord with the r absorption line -, while the other 
pair is seen at a strong redshift (500-900 km s^'). The ratio/? = F(f)/F(i) of the line fluxes is usu- 



ally used for density diagnostics [17]. However, in the presence of strong UV-radiation fields (such 
as in the vicinity of HDE 226868), the R ratio is systematically decreased due to photodepletion of 
the ls2s(^Si) level [13]. Therefore, an absolute density analysis would overestimate the density. 



One can, however, still infer the relative statement that the red-shifted pair of i- and f-lines (with 
the lower /?-ratio) stems from a much denser emitting plasma component. Both its velocity and its 
density indicate that this component can be identified with the focused wind. 

1.3 The dip-spectrum 

During the dips (Fig. |5|a), the soft X-ray spectrum is considerably absorbed. Defining the 
soft X-ray color as the ratio of the 0.5-1.5 keV and the 1.5-3 keV Chandra count rate, and the 
hard color as the ratio of 1.5-3 keV and 3-10 keV rate, a color-color diagram (Fig. ^) shows that 
neutral photoabsorption alone cannot explain the data. The nose-like shape of the data points can 
be explained by a partial covering model for the photoelectric absorption: 

(1 _/J .e-^/^r-M£) + .e-WH.2-M£) (1.3) 

We use a constant = 5.4x 10^^ cm^^ and a power-law continuum with photon index F = 1.73 
as found with RXTE-FCA for the non-dip spectrum [^. In this model, the strongest dips can be 
described by an additional absorber with Nu,2 > 4 x 10^^ cm^^ covering a fraction /c between 95 % 
and 98 % of the X-ray source. 

The high-resolution dip spectra (see, e.g.. Fig. |6|) show strong absorption lines of lower ionized 
ions, e.g., resonance Ka absorption lines of Si VI-XI between 6.8 and 7.1 A, which are not detected 
in the non-dip spectrum at all. This result is consistent with an origin of the dips in colder matter 



in the focused wind which covers the X-ray source. (More details will be presented in [ 10].) 
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Figure 5: (a): Qiandra-HETGS (0.5-10 keV) light curve (see also Fig. [l|). (b): Color-color diagram from 
(0.5-1.5keV):(1.5-3keV) and (1.5-3keV):(3-10keV) count rate ratios in 25 s intervals, explained by a grid 
of partial covering models (Eq. 1.3), assuming a power-law continuum (F = 1.73). The gray labels indicate 
the covering fraction /c in % and the column density A^h,2 of the additional absorber in 10^^ cm^^. 
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Figure 6: Chandra-HETGS Cyg X-1 specti-um of the first dip (^orb < 0.96, see Fig.il). While the non-dip 
specti'um (Fig. ||) showed silicon resonance Ka absorption Hnes of only Si XIV and XIII (Si Lya at 6.18 A 
and Si Hea at 6.65 A), the dip spectrum also presents the lower ionization stages Si VI-XI (green labels). 
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Table 1: Multi-satellite observation of CygnusX-1 on 2008 April 18/19. 



Satellite 



Start 



Stop 



Notes 



XMM 
INTEGRAL 



RXTE 
Suzaku 

Chandra 



{ 



{ 



2008-04-18 
2008-04-18 
2008-04-18 
2008-04-19 
2008-04-18 
2008-04-18 
2008-04-19 



12:07:27 2008-04-19, 04:25:46 

12:35:32 2008-04-19, 23:00:45 

13:21:36 2008-04-19,04:22:39 

17:37:20 2008-04-19, 20:04:47 

16:22:33 2008-04-19, 10:08:08 

18:27:19 2008-04-19,02:48:58 

14:57:27 2008-04-19,20:17:26 



C {obs. of the Cyg-region} 
1^ with SAA-interruptions 
with SAA-interruptions 



Note that, according to [|], 0orb(2OO8-O4-18, 15:53) = and (/)orb(2008-04-19, 20:00) = 0.20 , 



2. The multi-satellite campaign of 2008 April 

For a holistic study of CygX-1 and the focused wind of HDE 226868, further multi-satellite 
observations had to be performed. The strong variability due to dipping requires to observe Cyg X-1 
strictly simultaneously, if one wants to benefit from several different instruments in order to study 
the wind and the absorption dips. We originally proposed for a joint XMM+Chandra^ observation 
- with XMM's EPIC-pn camera for a high sensitivity for the spectral continuum shape through 
the broad iron line region, Chandra's HETGS for high-resolution wind spectroscopy, and XMM's 
RGS for a high-resolution spectrum at energies below the 0-edge. RXTE was available to provide 
the broadband continuum spectrum up to 250 keV. We were also able to obtain simultaneous INTE- 
GRAL time - extending the spectral coverage up to ~2 MeV -, as well as Suzaku - allowing for 
an independent confirmation of the spectral continuum, see Table |l|. We were furthermore awarded 
a few ks of simultaneous Swift TOO time. AGILE did not detect the source in the 30 MeV-50 GeV 
band. Cyg X-1 has also been monitored in the radio by the VLA during this observation. 

As expected, the observation was again dominated by dipping events. Most of them are re- 
stricted to the soft X-ray band (Chandra, XMM), but a few are also detected above 4keV by the 
PCA on RXTE or even above 20keV by IBIS on INTEGRAL. A color-color diagram reveals the 
same nose-like shape as was found for the dips in 2003 (Fig. ^), indicating the same possible 
origin of the dips in partial covering of the X-ray source [ |l6| , |lO| ]. 

Due to the fast decrease of the photoelectric absorption cross section Obf above the ionization 
threshold. X-ray photons with energies well above the iron edge at 7. 1 keV can hardly be absorbed 
by the wind. They can, however, still be scattered out of the line of sight. Figure |7| shows such an 
X-ray scattering trough lasting more than 8 hours: Between April 18, ~16:45 and April 19, ~1:15, 
the 20-40 keV count rate measured with INTEGRAL -IBIS, as well as the 12-60keV count rate 
measured with Suzaku-PIN, decrease until April 18, 20:15 and rise again afterwards. A continuous 
model joining two linear functions (a descending one followed by an ascending one) can describe 
both light curves simultaneously (with a scaling factor of 0.23 between the PIN and IBIS count 
rates) remarkably well (Z^^^ = 2). If the 30 % reduction between the highest and the lowest count 
rate of this fit (corresponding to non-dip and the strongest dip) is due to Compton scattering, the 
optical depth T = 0.36 corresponds to Nu = 6 x lO^^cm^^. 

^ Note that scheduling an simultaneous observation with orbital phase constraints is not easy for these two satellites... 
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Figure 7: Light curves obtained with INTEGRAL-lBlS in the 20-40 keV band (top black data, left y axis) 
and Suzaku-PYN in the 12-60keV band (lower blue data, right y-axis) with a time resolution of 200 and 96 s, 
revealing an X-ray scattering trough shortly after superior conjunction, which lasts for 8.5 hours. During 
this time, the light curves can be described by two continuously jointed linear functions (red lines). 




Energy [keV] 

Figure 8: The 0.5keV-2MeV broadband spectrum of CygX-1, observed simultaneously with XMM, 
Suzaku, RXTE and INTEGRAL in 2008 April. The data below 300 ke V can be described consistently by 
a broken power-law with exponential cutoff at high energies (fifoid ~ 170keV). ISIS' (model independent) 
flux-corrected spectra (relying on the response matrices) seem to indicate a non-thermal tail in the SPI data. 



The broadband spectrum between 0.5 keV and 2 MeV, assembled from various instruments 
(Fig. |8|), reveals that CygX-1 was again observed in the hard state. Below 300 keV, it can be 
described as a simple elbow-shaped broken power-law with high energy cutoff. RXTE-PCA's 
photon indices below and above f'break = 9.8keV, namely Fi = 1.7 and F2 = 1.4, are amongst the 
hardest ones found in the long term monitoring of CygX-1 from 1999 to 2004 [1^ 18]. Preliminary 
INTEGRAL-SPl data seem to indicate a non-thermal tail in the spectrum above 400 keV, but further 
analysis through joint broadband fits with Comptonization models is still ongoing. 
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